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ABSTRACT 

Aims. Dust-obscured star-formation becomes much more important with increasing intensity, and increasing redshift. We aim to 
reveal cosmic star-formation history obscured by dust using deep infrared observation with the AKARI. 

Methods. We construct restframe 8/im, 12/im, and total infrared (TIR) luminosity functions (LFs) at 0.15 < z < 2.2 using 4128 
infrared sources in the AKARI NEP-Deep field. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, 
and 24/im) by the AKARI satellite allows us to estimate restframe 8/im and 12/^m luminosities without using a large extrapolation 
based on a SED fit, which was the largest uncertainty in previous work. 

Results. We have found that all 8/im (0.38 < z < 2.2), 12/im (0.15 < z < 1.16), and TIR LFs (0.2 < z < 1.6), show a continuous 
and strong evolution toward higher redshift. In terms of cosmic infrared luminosity density (S1/h), which was obtained by integrating 
analytic fits to the LFs, we found a good agreement with previous work at z < 1.2. We found the evolves as oc (1 + a)'* ''*^ *'. 
When we separate contributions to Qir by LIRGs and ULIRGs, we found more IR luminous sources are increasingly more important 
at higher redshift. We found that the ULIRG (LIRG) contribution increases by a factor of 10 (1.8) from z=0.35 to 2=1.4. 

Key words, galaxies: evolution, galaxies:interactions, galaxies :starburst, galaxies :peculiar, galaxies:formation 

1 . Introduction iBell et alJ (l2005h estimate that IR luminosity density is 7 

times higher than the UV lu minosity density at z^O.7 than lo- 

Studies of the extragalactic background suggest at least half ^^jj Xakeuchi. BuaL^Burgarella (2005) reported that UV-to- 

the luminous energy generated by stai's has been reprocessed ir i™osity density ratio, evolves from 3.75 

intothe infrared (IR) by dust (Lagache et al., 1999; Puget et al., ^^^q^ ^5 ^ q ^-^^^ ^ ^^^f^j \^^^^^^^^ sample 

[19961 i Franceschim^Rodighiero. & Vaccan, 2008), suggesting selection effect, and that 70% of star formation activity is ob- 

that dust-obscured star formation was much more important at ^^^^^^ ^^^^ q 5< ^^12. Both works highlight the im- 

higher redshifts than today. portance of probing cosmic star formation activity at high red- 

~ shift in the infrared bands. Several works found that most ex- 

This research is based on the observations with AKARI, a JAXA . . ^ ■ /ct^ ^ ■ u-u • -i- 

... . . . treme star-forming (SF) galaxies, which are increasingly impor- 

proect with the participation ot ESA. ,. , , , •, , , 

. . ^ „.j.oi, T- , uu- .ji, tant at higher redshifts, are also more heavily obscured by dust 

Based on data collected at Subaru Telescope, which is operated by itt — . . " 71 1 -„„ ,1 1,, ,,. 71 1-„„ - i r=r — -1 rrfi_, 

the National Astronomical Observatory of Japan. ( |Hopkms et aU |200ll| Sulhvan et al-lM | Buat et al. L|2007D. 

JSPS SPD fellow 
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Despite the value of infrared observations, studies of 
infrared galaxies by the IRAS and the ISO were re- 
stricted to bright so urces due to the limited sensitiv- 

ities jSa unders et^ 1 19901: iR owan-Robinson et al., 199J 
Flore s et al., 1999'; 'Serieant et al., 2004; .Takeuchi et al.. .200" 
Take uchi. Yoshikawa, & Ishii, 2003), until the recent launch of 
the Spitzer and the AKARI satellites. Their enormous improved 
sen sitivities have revolutio nized the field. For example: 

iLe Floc'h et al.l (l2005l) analyzed the evolution of the total 
and 15/im IR luminosity functions (LFs) at < z < 1 based 
on the the Spitzer MIPS 24^m data (> 83//Jy and R < 24) in 
the CDF-S, and found a positive evolution in both luminosity 
and density, suggesting increasing importance of the LIRG and 
UL IRG populations at higher re dshifts . 

iPerez-Gonzalez et~aL (l2005h used MIPS 24/im observations 
of the CDF-S and HDF-N (> 83/iJy) to find that that L* steadily 
increases by an order of magnitude to z ^ 2, suggesting that the 
luminosity evolution is stronger than the density evolution. The 
nx iR scales as (l4- z )^ °^° -^ from z=Q to 0.8. 

iBabbedg e et"an (l2006h constructed LFs at 3.6, 4.5, 5.8, 8 
and 24/im over < z < 2 using the data from the Spitzer 
Wide-area Infrared Extragalactic (SWIRE) Survey in a 6.5 deg^ 
(S24tj.m > 230/iJy). They found a clear luminosity evolu- 
tion in all the bands, but the evolution is more pronounced at 
longer wavelength; extrapolating from 24/im, they inferred that 
^TiR oc(Ih-z)^'^. They constructed separate LFs for three dif- 
ferent galaxy SED (spectral energy distribution) types and Type 
1 AGN, finding that starburst and late-type galaxies showed 
stronger evolution. Comparison of 3.6 and4.5/imLFs with semi- 
analytic and spectrophotometric models suggested that the IMF 
is skewed towards higher mass star formation in more intense 
sta rbursts. 

ICaputi et al.l (l2007l) estimated restframe 8/im LFs of galaxies 
over O.OSdeg^ in the GOODS fields based on Spitzer 24/im (> 
80/iJy) at z=l and 2. They found a continuous and strong posi- 
tive luminosity evolution from z=0 to z=l, and to z=2. However, 
they also found that the number density of star-forming galaxies 
with j/L^^™ > 10^°-^Lo (AGNs are excluded.) increases by a 
factor of 20 from z=Q to 1, but decreases by half from z=l to 2 
ma inly due to the decreas e of LIRGs. 

iMagnelh et all (l2009l) investigated restframe 15/im, 35/im 
and total infrared (TIR) LFs using deep 70/im observations 
(~300 /iJy) in the Spitzer GOODS and FIDEL (Far Infrared 
Deep Extragalactic Legacy Survey) fields (0.22 deg^ in total) 
at z < 1.3. They stacked 70/im flux at the positions of 24/im 
sources when sources are not detected in 70/im. They found no 
change in the shape of the LFs, but found a pure luminosity evo- 
lution proportional to (1h-z)^ and that LIRGs and ULIRGs 
have increased by a factor of 40 and 100 in number density by 
z~l. 

Also, see lDai et al.] (12009 ') for 3.6-8.0 /im LFs based on the 
IRAC photometry in the NOAO Deep Wide-Field Survey Bootes 
field. 

However, most of the Spitzer work relied on a large 
extrapolation from 24/im flux to estimate the 8, 12/im or 
TIR luminosity. Consequently, Spitzer results heavil y de- 
pended on the assumed IR SED li brary (Dale & Helou, I2OOI 
iLagache. Dole. & Pugeli 120031: ICharv & Elbazi 120011) . Indeed 
many authors pointed out that the largest uncertainty in these 
previous IR LFs came from SED models, especially when one 
computes TIR luminosity solely from observed 24/im flux (e.g., 
see Fig. 5 of Caputi et al., 2007). 

AKARI, the first Japanese IR dedicated satellite, has con- 
tinuous filter coverage across the mid-IR wavelengths, thus, al- 




Fig. 1. P h otome t ric redshift e stimate s with Le P hare 
(Ilb ertetaLL I2OO6I: lArnouts et all IIpOTL lllbert et al.L I2OO9I) 
for spectroscopically observed galaxies with Keck/DEIMOS 
(Takagi et al. in prep.). Red squares show objects where AGN 
templates were better fit. Errors of the photoz is ^^^=0.036 for 

z < 0.8, but becomes worse at z > 0.8 to be ^^^=0.10 due 
mainly to the relatively shallow near-IR data. 



lows us to estimate MIR (mid-infrared)-luminosity without us- 
ing a large fc-correction based on the SED models, eliminating 
the largest uncertainty in previous work. By taking advantage of 
this, we present the restframe 8, 12/im and TIR LFs using the 
AKARI NEP-Deep data in this work. 

Restframe 8/im luminosity in particular is of primary rele- 
vance for star-forming galaxies, as it includes polycyclic aro- 
matic hydrocarbon (PAH) emiss ion. PAH molecules charac - 
terize star-forming regions (Des ertTBoulanger. & Pugelll990l) . 
and the associated emission lines between 3.3 and 17 /im dom- 
inate the SED of star-forming galaxies with a main bump lo- 
cated around 7.7/im. Restframe 8/im luminosities have been 
confirmed to be go od indicators of knots of star formation 
dCalzetti et al.L l2005h and of the overall s tar formation activity 
of star forming galaxies dWu et all l2005h . At z=0.375, 0.875, 
1.25 and 2, the restframe 8/im is covered by the AKARI S'll, 
L15, L18W and L24 filters. We present the restframe 8/im LFs 
at these redshifts at Sectio HTT] 

Restframe 12/im lumi nosity functions have a lso been 
studied extensively (R ush. Malkan. & Spinogliol 119931 ; 
Perez-Gonzalez et al., 2005]). At z=0.25, 0.5 and 1, the 
restframe 12/im is covered by the AKARI L15, L18W and L24 
filters. We present the restframe 12/im LFs at these redshifts in 
Sectior]33] 

We also estimate TIR LFs through the SED fit using all 
the mid-IR bands of the AKARI. The results are presented in 
Section[33] 

Unless otherwise stated, we adopt a cosmology with 
(/i, rim, ^a) = (0.7, 0.3, 0.7) dKomatsu et al.Ll2008) . 

2. Data & Analysis 

2. 1. Multi-wavelength data in ttie AKARI NEP Deep field 

AKARI, the Japanese infrared satellite (Murak ami et al.Ll2007h . 
performed deep imaging in the North Ecliptic Region (NEP) 
from 2-24/im, with 1 4 pointings in each field ov er 0.4 
deg2 dMatsuharaetall l2006l 120071: IWada etall l2008h . Due 
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Fig. 2. Photometric redshift distribution. 
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Fig. 3. 8/im luminosity distributions of samples used to compute 
restframe 8/im LFs. From low redshift, 533, 466, 236 and 59 
galaxies are in each redshift bin. 



to the solar synchronous orbit of the AKARI, the NEP 
is the only AKARI field with very deep imaging at these 
wavelengths. The 5 a sensitivity in the AKARI IR filters 
(N2, m, N4, S7, S9W, Sll, L15, L18W and L24) are 14.2 , 
11.0, 8.0, 48, 58, 71, 117, 121 and 275 uJv (IWada et al.L l2008h . 
These filters provide us with a unique continuous wavelength 
coverage at 2-24/im, where there is a gap between the Spitzer 
IRAC and MIPS, and th e ISO LW2 and LV^3. Please consult 
IWada et al.l (l2007l l2008l) : IPearson et all (l2009albl) for data veri- 
fication and completeness estimate at these fluxes. The PSF sizes 
are 4.4, 5.1, and 5.4" in 2 - 4, 7 - 11, 15 - 24/im bands. The 
depths of near-IR bands are limited by source confusion, but 
those of mid-IR bands are by sky noise. 



In analyzing these observations, we first combined the three 
images of the MIR channels, i.e. MIR-S(57, S'QVF, and 5*11) 
and MIR-L(L15, L18W and L24), in order to obtain two high- 
quality images. In the resulting MIR-S and MIR-L images, the 
residual sky has been reduced significantly, which helps to ob- 
tain more reliable source catalogues. For both the MIR-S and 
MIR-L channels, we use SExtractor for the combined images to 
determine initi al source posi t ions. 

We follow iTakagi et al.l (l2007h procedures for photometry 
and band-merging of IRC sources. But this time, to maximize 
the number of MIR sources, we made two IRC band-merged 
catalogues based on the combined MIR-S and MIR-L images, 
and then concatenated these catalogues, eliminating duplicates. 

In the band-merging process, the source centroid in each IRC 
image has been determined, starting from the source position in 
the combined images as the initial guess. If the centroid deter- 
mined in this way is shifted from the original position by > 3", 
we reject such a source as the counterpart. We note that this 
band-merging method is used only for IRC bands. 

We compared raw number counts with previous work based 
on the same data but with different sourc e extraction methods 
(IWada et al.L 120081; iPearson et al.L l2009allbl) and found a good 
agreement. 

A subregion of the NEP-Deep field wa s observed in the 
BVRi'z'-h mds with the Subaru telescope ("Imai e t aTl 120071: 
IWada et all 12008) . reaching limiting magnitudes of zab =26 
in one field of view of the Suprime-Cam. We restrict our analy- 
sis to the data in this Suprime-Cam field (0.25 deg^), where we 
have enough UV-opical-NIR coverage to estimate good photo- 
metric redshifts. The u'-band photometry in this area is provided 
by the CFHT (Serjeant et al. in prep.). The same field was also 
observed with the KPN 02m/FLAMINGO s in J and Ks to the 
depth of Ksvega < 20 (llmaiet al.Ll2007h . GALEX covered the 
entire field to depths of FUV < 25 and NUV < 25 (Malkan et 
al. in prep.). 

In the Suprime-Cam field of the AKARI NEP-Deep field, 
there are a total of 4128 infrared sources down to ~'100 /iJy in 
the L18W filter All magnitudes are given in AB system in this 
paper 

For the optical identification of MIR sources, we adopt th e 
UkeHhood ratio (LR) method ( Suflierland & Saundersl 1 19921) . 
For the probability distribution functions of magnitude and an- 
gular separation based on correct optical counterparts (and for 
this purpose only), we use a subset of IRC sources, which are 
detected in all IRC bands. For this subset of 1100 all-band- 
detected sources, the optical counterparts are all visually in- 
spected and ambiguous cases are excluded. There are multiple 
optical counterparts for 35% of MIR sources within < 3". If we 
adopted the nearest neighbor approach for the optical identifica- 
tion, the optical counterparts differs from that of the LR method 
for 20% of the sources with multiple optical counterparts. Thus, 
in total we estimate that less than 15% of MIR sources suffer 
from serious problems of optical identification. 

2.2. Photometric redshift estimation 

For these infrared sources, we have computed photomet- 
ric redshift using a pubhcly av ailable code, LeP h are P] 
(lllbert et all llooa lArnouts et all HMk lllbert et all I2009Ik 
The input magnitudes are FUV, NUV(GALEX),u(CFUT), 
B, V, R, i', z'(Subaru), J, and if (KPN02m). We summarize the 
filters used in Table [T] 

' http://www.cfht.hawaii.edu/~arnouts/lephare.html 
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Table 1. Summary of filters used. 


Estimate 


Redshift 


Filter 


Photo z 


0.15<z<2.2 


FUV, NUV,u, B, V, R, i',z, J, and K 


8/im LF 
8/im LF 
8/im LF 
8/im LF 


0.38<z<0.58 
0.65<z<0.90 
l.l<z<1.4 
1.8<z<2.2 


SU (U/im) 
L15 (15/tm) 
L18W (18/im) 
L24 (24/im) 


12/tm LF 
12/im LF 
12/im LF 


0.15<z<0.35 
0.38<z<0.62 
0.84<z<1.16 


L15 (15/im) 
L18W (18/im) 
L24 (24/im) 


TIRLF 
TIR LF 
TIR LF 
TIR LF 


0.2<z<0.5 
0.5<z<0.8 
0.8<z<1.2 
1.2<z<1.6 


S7, S9W, Sll, L15, L18W and L24 
S7, S9W, Sll, L15, L18W and L24 
S7, S9W, Sll, L15, L18W and L24 
S7, S9W, Sll, L15, L18W and i24 



redshift of the redshift bin. More precisely, 

Zmax = niin(zniax of thc bin, Zniax from the flux limit) 



(1) 



Among various templates and fitting parameters we tried, 
we found the best resul ts were obtained with the fo l lowing : we 
used modified CWW jColeman. Wu. & Weedmanl \mat) and 
QSO templates. These CWW templat es are interpolated an d ad- 
justed to better match VVDS spectra (Arnouts et al.L l2007h . We 
included strong emission lines in computing colors. We used 
the Calzet ti extinction law. M ore details in training LePhare 
is given in lllbert et aP (l2006l) . 

The resulting photometric redshift estimates agree reason- 
ably well with 293 galaxies (R < 24) with spectroscopic red- 
shifts taken with Keck/DEIMOS in the NEP field (Takagi et al. 
in prep.). The measured errors on the photo-z are ^^=0.036 for 

z < 0.8 and 4^=0.10 for z > 0.8. The ^ becomes signifi- 

— 1+2 1+Z ^ 

cantly larger at z > 0.8, where we suffer from relative shallow- 
ness of our near-IR data. The rate of catastrophic failures is 4% 
(^^^ >0.2) among the spectroscopic sample. 

In FiglT] we compare spectroscopic redshifts from 
Keck/DEIMOS (Takagi et al.) and our photometric red- 
shift estimation. Those SEDs which are better fit with a QSO 
template are shown as red triangles. We remove those red 
triangle objects (~2% of the sample) from the LFs presented 
below. We caution that this can only remove extreme type-1 
AGNs, and thus, fainter, type-2 AGN that could be removed by 
X-rays or optical spectroscopy still remain in the sample. 

Fig|2]shows the distribution of photometric redshift. The dis- 
tribution has several peaks, w hich corresponds to galaxy clusters 
in the field dGotoet all 12008 ). We have 12% of sources that do 
not have a good SED fit to obtain a reliable photometric redshift 
estimation. We apply this photo-z com pleteness correction to th e 
LFs we obtain. Readers are referred to lNegrello et at aL I (l2009l) . 
who estimated photometric redshifts using only the AKARI fil- 
ters to obtain 10% accuracy. 



2.3. The Wrn^^ method 



We compute LFs using the l/14iax method (ISchmidtlll968h . The 
advantage of the 1/Knax method is that it allows us to compute 
a LF directly from data, with no parameter dependence or an 
assumed model. A drawback is that it assumes a homogeneous 
galaxy distribution, and is thus vulnerable to local over-/under- 
densities (Takeuchi. Yoshikawa. & IshiL .2000). 

A comoving volume associated with any source of a given 
luminosity is defined as V,nax — ^z„,ax ~ ^z„,i„, where Zmin 
is the lower limit of the redshift bin and z„iax is the maximum 
redshift at which the object could be seen given the flux limit of 
the survey, with a maximum value corresponding to the upper 



We used the SED templates d Lagache. Dole. & Pugei.l2003h for 
fc-corrections to obtain the maximum observable redshift from 
the flux limit. 

For each luminosity bin then, the LF is derived as 



1 



-Wi, 



(2) 



where l^iax is a comoving volume over which the ith galaxy 
could be observed, AL is the size of the luminosity bin (0.2 dex), 
and Wi is the completeness correction facto r of the ith galaxy . 
We use completeness correction measured by Wada et al. ("2008) 
for 1 1 and 24 /tm and Pearson et al. (2009a b) for 15 and 18/im. 
This correction is 25% at maximum, since we only use the sam- 
ple where the completeness is greater than 80%. 

2.4. Monte Carlo simulation 

Uncertainties of the LF values stem from various factors such 
as fluctuations in the number of sources in each luminosity bin, 
the photometric redshift uncertainties, the /c-correction uncer- 
tainties, and the flux errors. To compute these errors we per- 
formed Monte Carlo simulations by creating 1000 simulated cat- 
alogs, where each catalog contains the same number of sources, 
but we assign each source a new redshift following a Gaussian 
distribution centered at the photometric redshift with the mea- 
sured dispersion of Az/(1 + z) =0.036 for z < 0.8 and 
Az/(1 + z) =0.10 for z > 0.8 (Fig[B. The flux of each source 
is also allowed to vary according to the measured flux error fol- 
lowing a Gaussian distribution. For 8/im and 12/im LFs, we can 
ignore the errors due to the /c-correction thanks to the AKARI 
MIR filter coverage. For TIR LFs, we have added 0.05 dex of 
error for uncertainty in the SED fitting following the discussion 
in Magnelli et al. (2009). We did not consider the uncertainty 
on the cosmic variance here since the AKARI NEP field cov- 
ers a large vol ume and has comparab le number counts to other 
general fields dlmai et all I2007l l2008l). Each redshift bin we use 
covers ^10^ Mpc^ of volume. See'Mats uhara et aTl (1200 6) for 
more discussion on the cosmic variance in the NEP field. These 
estimated errors are added to the Poisson errors in each LF bin 
in quadrature. 



3. Results 

3.1. Slim LF 

Monochromatic 8/im luminosity (Ls um) is known t o cor- 
relate well with the TIR luminosity iBabbedpe et all 120061: 
iHuang et al.Ll2007 ). especially for star-forming galaxies because 
the rest-frame 8/im flux are dominated by prominent PAH fea- 
tures such as at 6.2, 7.7 and 8.6 /im. 

Since the AKARI has continuous coverage in the mid-IR 
wavelength range, the restframe 8/im luminosity can be obtained 
without a large uncertainty in /c-correction at a corresponding 
redshift and filter. For example, at z=0.375, restframe 8/im is 
redshifted into 511 filter Similai'ly, Lib, L18W and L24 cover 
restframe 8/im at z=0.875, 1.25 and 2. This continuous filter 
coverage is an advantage to AKARI data. Often SED models are 
used to extrapolate from Spitzer 24/im flux in previous work. 
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producing a source of the largest uncertainty. We summarise fil- 
ters used in Table [T] 

To obtain restframe 8/im LF, we applied a flux limit 
of F(S11)<70.9, F(L15)<117, F(L18W)<121.4, and 
F(L24)<275.8 ^iJy at z=0.38-0.58, z=0.65-0.90, z=l. 1-1.4 
and z= 1.8-2.2, res pectively. These are the Sa limits measured in 
IWada et all 12008). We exclude those galaxies whose SEDs are 
better fit with QSO templates (©. 

We use the completeness cu rve presented in IWada et aP 
(120081) and iPearson et al.) (l2009allbl) to correct for the incom- 
pleteness of the detection. However, this correction is 25% at 
maximum since the sample is 80% complete at the Sa limit. Our 
main conclusions are not affected by this incompleteness correc- 
tion. To compensate for the increasing uncertainty in increasing 
z, we use redshift binsize of 0.38< z <0.58, 0.65< z <0.90, 
1.1< z <1.4, and 1.8< z <2.2. We show the Ls^m distribution 
in each redshift range in Fig|3] Within each redshift bin, we use 
1/Vmax method to compensate for the flux limit in each filter. 

We show the computed restframe 8/im LF in Fig|4] Arrows 
show the 8/^m luminosity corresponding to the flux limit at the 
central redshift in each redshift bin. Errorbars on each point are 
based on the Monte Carlo simulation (i) l2.3l l. 

For a comparison, as the green dot-dashed line, we also 
sh ow the 8Mm LF of star-forming galaxies at < z < 0.3 
bv lHuan^ et al.l (l2007h . using the 1/Vmax method applied to the 
IRAC 8/xm GTO data. Compared to the local LF, our 8/xm LFs 
show strong evolution in luminosity. In the range of 0.48 < z < 
2, -Lg^nj evolves as oc (1 + z)^-^**'-^. Detailed comparison with 
the hterature will be presented in ^ 

3.2. Bolometric IR luminosity density based on tine 8^m 
LF 

Constraining the star formation history of galaxies as a func- 
tion of redshift is a key to understanding galaxy formation in 
the Universe. One of the primary purposes in computing IR 
LFs is to estimate the IR luminosity density, which in turn is a 
good estimator of the dust hidden cosmic star formation density 
dKennicuttl [19981) . Since dust obscuration is more important for 
more actively star forming galaxies at higher redshift, and such 
star formation cannot be observed in UV light, it is important to 
obtain IR-based estimate in order to fully understand the cosmic 
star formation history of the Universe. 

We estimate the total infrared luminosity density by integrat- 
ing the LF weighted by the luminosity. First, we need to convert 
Lsfim to the bolometric infrared luminosity. The bolometric IR 
luminosity of a galaxy is produced by the thermal emission of 
its interstellar matter. In star-forming galaxies, the UV radiation 
produced by young stars heats the interstellar dust, and the repro- 
cessed light is emitted in the IR. For this reason, in star-forming 
galaxies, the bolometric IR luminosity is a good estimator of the 
curren t SFR (star formation rate) of the galaxy. iBavouzet et al.l 
(120081) showed a strong correlation between Lg^m and total in- 
frared luminosity (Lti r) for 372 lo cal star-forming galaxies. 
The conversion given bv lBavouzet et al.. (.2008.) is: 



^TIR 



377.9 X (z.L.)0fJg^,„(±37%) 



(3) 



ICaputi et al.l (l2007t) further constrained the sample to lumi- 
nous, high S/N galaxies (i^i^^™ > 10^°Lq and S/N> 3 in all 
MIPS bands) in order to better match their sample, and derived 
the following equation. 




X 1 .8<z<2.2 
□ 1 .1<z<1 .4 
A0.65<z<0.90 
O0.38<z<0.58 



9.5 10.0 10.5 1 1.0 1 1.5 12.0 
log vLj8iJ,m') 



Fig. 4. Restframe 8/im LFs based on the AKARI NEP-Deep 
field. The blue diamons, purple triangles, red squares, and or- 
ange crosses show the 8/im LFs at 0.38 < z < 0.58, 0.65 < 
z < 0.90,1.1 < z < 1.4, and 1.8 < z < 2.2, respectively. 
AKARJ's MIR filters can observe restframe 8/im at these red- 
shifts in a corresponding filter Errorbars are from the Monte 
Caro simulations (i ]2.4! ). The dotted lines show analytical fits 
with a double-power law. Vertical arrows show the 8/im lumi- 
nosity corresponding to the flux limit at th e centr al redshift in 
each redshift bin. Overplotted are Babbed ge et al.l (j2006) in the 
pink dash-dotted hnes. jCaputi et al.. (200% in the cyan dash- 
dotted lines, and iHuang et al.l ( 2007h in the green dash-dotted 
lines. AGNs are excluded from the sample ( ^2.2i . 



Ltir. = 1-91 x {vL,)l-^!,^^^{±bb%) (4) 

Since ours is also a sample of bright galaxies, we use this 
equation to convert L^^m to Ltir. Because the conversion is 
based on local star-forming gal axies, it i s a co ncern if it holds 
at higher redshift or not. iBavouzet et al.l (|2P08') checked this by 
stacking 24/im sources at 1.3 < z < 2.3 in the GOODS fields 
to find the stacked sources are consistent with the local rela- 
tion. They concluded that equation (|3]l is v alid to link ig/jm 
and Ltir at 1.3 < z < 2.3. Takagi et al.l (12010) also show 
that local ^7.7^^™ vs Lti r relation h olds true for IR galaxies 
at Z'-l (see their Fig. 101. IPope et all (i2008) showed that z --2 
sub-millimeter galaxies lie on the relation between Ltir and 
LpAH,7.7 that has been established for local starburst galaxies. 
570/5*24 ratios of 70/im sources in Papovich et al. (2007) are 
also consistent with local SED templates. These results suggest 
it is reasonable to use equation (|4|i for our sample. 

The conversion, however, has been the largest source of er- 
ror in estimating Ltir from Lgp,„. Bavouzet et al. (2008) them- 
selves quote 37% of uncertainty, and that Caputi et al. (2007) re- 
port 55% of dispersion around the relation. It should be kept in 
mind that the restframe 8/im is sensitive to the star-formation 
activity, but at the same time, it is where the SED models have 
strongest discrepancies due to the complicated PAH emission 
lines. A detailed c omparison of diffe rent conversions is pre- 
sented in Fig. 12 of ICaputi et al.l (|2007|) . who reported factor of 
^^5 of differences among various models. 
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Then the 8^m LF is weighted by the Ltir and integrated 
to obtain TIR density. For integration, we first fit an ana- 
lytical function to the LFs. In the literature, IR LFs were 
fit better by a doubl e-power law d Babbedge et al., 2006) or 
a double-expon ential (Saunders et al.(, ,1990|: .Pozzi et al.> ,2004t 
iTakeuchi et all [2006; Le Floc'h et all l2005h than a Schechter 
function, which declines too suddenlly at the high luminosity, 
underestimating the number of bright galaxie s. In this work , 
we fit the 8 /im LFs using a double-power law ( jBabbedge et al.L 
I2006h as shown below. 



^{L)dL/L* = 



^[L)dL/L* ^ 



dL/L*, {L < L* 



1-/3 



J7 ) dL/L* 



(L > L*) 



(5) 



(6) 



G 

a 



0.0 



This work 

Le Floth et al. 



.00 



0.10. 







0.01 



0.5 1.0 1.5 2.0 

z 



First, the double-power law is fitted to the lowest redshift LF at 
0.38< z <0.58 to determine the normalization ($*) and slopes 
(a, P). For higher redshifts we do not have enough statistics to si- 
multaneously fit 4 parameters ($*, L*, a, and (3). Therefore, we 
fixed the slopes and normalization at the local values and varied 
only L* at for the higher-redshift LFs. Fixing the faint-end slope 
is a c ommon procedure with t he depth of current IR satellite sur- 
veys dBabbedge et all 120061: ICaputi et all l2007h . The stronger 
evolution in luminosity than in density found by previous work 
(Perez-Gonzalez et al., 2005; Le Floc'h et al., 2005) also justi- 
fies this parametrization. Best fit parameters are presented in 
Table 12] Once the best-fit parameters are found, we integrate the 
double power law outside the luminosity range in which we have 
data to obtain estimate of the total infrared luminosity density. 

The resulting total luminosity density (fi/ij) is shown in 
Fig|5]as a function of redshift. Errors are estimated by varying 
the fit within Icr of uncertainty in LFs, then errors in conversion 
from L^^m to Ltir are added. The latter is by far the larger 
source of uncertainty. Simply switching from equation (O (or- 
ange dashed line) to (|4|i (red solid line) produc es a ^50% dif- 
ferenc e. Cyan dashed lines show results from iLe Floc'h et al.l 
(I2005h for a co mparision. The lowes t redshift point was cor- 
rected following lMagnelh et aP (l2009l) . 

We also show the evolution of monochromatic 8/im lumi- 
nosity (Lsfim), which is obtained by integrating the fits, but 
without converting to Ltir in Fig|6l The ilsfj.m evolves as 
oc (1 + 

The SFR and Ltir. are related by the following equation 
for a Salpeter IMF, (p (m) cx m"^-^^ between 0.1 - IOOMq 
d^nicutt, 1998). 

SFR{MQyr-') - 1.72 x 10' Ltir{Lq) (7) 

The right ticks of Fig|5] shows the star formation density 
scale, converted from fliR using the above equation. 

In Fig|5] iliR monotonically increases toward higher z. 
Compared with z=Q, fliR is ~10 times larger at z=l. The evolu- 
tion between z=0.5 and 2;= 1 .2 is a little flatter, but this is perhaps 
due to a more irregular shape of LFs at 0.65 < z <0.90, and thus, 
we do not consider it significant. The results obtained here agree 
with previous work (e.g., LeJilQclh et al., 2005) within the er- 
rors. We compare the results with previous work in more detail 
ing] 



Fig. 5. Evolution of TIR luminosity density computed by inte- 
grating the 8/im LFs in Fig|4]The red solid lines use the con- 
version in equation (|4|. T he orange da shed lines use equation 
(O. Results from lLe Floc'h et alj(l2005l) are shown with the cyan 
dotted lines. 
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Fig. 6. Evolution of 8/im IR luminosity density computed by 
integrating the 8/im LFs in Fig|4] The lowest redshift point is 
from lHuang et all (120071) . 



3.3. 12fimLF 

In this subsection we estimate restframe 12/im LFs based 
on the AKARI NEP-Deep data. 12/im luminosity {Li2^m) 
has been well studied through ISO an d IRAS, and known to 
correlate closely with TIR luminosity (ISpinoglio et al.L 119951; 
iPerez-Gonzalez et al.Ll2005h . 

As was the case for the 8/im LF, it is advantageous that 
AKARI's continuous filters in the mid-IR allow us to estimate 
restframe 12/im luminosity without much extrapolation based 
on SED models. 
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Table 2. Best fit parameters for 8,12/im and TIR LFs 



Redshift 


A 




$*(Mpc"^dex 




a 




0.38<z<0.58 


8/im 


(2.2+l'l) X 10'" 


(2.11^,:^) X 10- 


3 


1 


7C+U.01 

—0.01 




0.65<z<0.90 


8/im 


(2.8+o'i) X 10'" 


2.1 X IQ-^ 




1 


75 


3.5 


l.l<z<1.4 


8/im 


(3.3+EJ-^) X lO'o 


2.1 X 10"^ 




1 


75 


3.5 


1.8<z<2.2 


8/im 


(8.2ii-2) X 10'° 


2.1 X 10"^ 




1 


75 


3.5 


0.15<z<0.35 


12/im 


(6.8+^-i) X 10« 


(4.21)^:'^) X 10- 




1 


^"-0.02 


7 Q+U.4 


0.38<z<0.62 


llfim 


(UJtti) X 10» 


4.2x10"^ 




1 


20 


2.9 


0.84<z<1.16 


12/im 


(14l^) X 10^ 


4.2x10-=^ 




1 


20 


2.9 


0.2<z<0.5 


Total 


(1.2+^,^) X 10" 


(5.6+02) X 10- 


4 


1 


o+U.i 
S-0.4 


3 0+'-" 


0.5<z<0.8 


Total 


(2.4j:}:*) X 10" 


5.6x10"^ 




1 


8 


3.0 


0.8<z<1.2 


Total 


(3.91:2:2) X 10" 


5.6x10"'' 




1 


8 


3.0 


1.2<z<1.6 


Total 


(14+^) X 10" 


5.6x10"" 




1 


8 


3.0 




Fig. 7. 12/im luminosity distributions of samples used to com- 
pute restframe 12/im LFs. From low redshift, 335, 573, and 213 
galaxies are in each redshift bin. 



Targeted redshifts are z=0.25, 0.5 and 1 where Lib, L18W 
and i24 filters cover the restframe 12/im, respectively. We sum- 
marise the filters used in Table [1] Methodology is the same as 
for the 8/im LF; we used the sample to the Scr limit, corrected 
for the completeness, then used the 1/Fmax method to com- 
pute LF in each redshift bin. The histogram of ii2pm distri- 
bution is presented in Fig|7] The resulting 12/xm LF is shown 
in FiglH] Compared with Rush. Malkan. & Spinogliol (jl993)'s 
z=Q LF based on IRAS Faint Source Catalog, the 12/im LFs 
show steady evolution with increasing redshift. In the range of 
0.25 < z < 1, i^sA^m evolves as oc (1 + z)^-^^°-'^. 

3.4. Bolometric IR luminosity density based on tine 12^im 
LF 

12/Lim is one of the most frequently used monochromatic fluxes 
to estimate Ltir. The total infrared l uminosity is compute d 
from the Li2am using the co nversion in ICharv&ElbaldlooTl) : 



iPerez-Gonzalez et 



iing tnecoi 
aiT (l2005h . 



Fig. 8. Restframe 12/im LFs based on the AKARI NEP-Deep 
field. The blue diamonds, purple triangles, and red squares show 
the 12/im LFs at 0.15 < z < 0.35,0.38 < z < 0.62, and 
0.84 < z < 1.16, respectively. Vertical arrows show the 12/im 
luminosity corresponding to the flux li mit at the central red- 
shift in each redshift bin. Overplotted are iPerez-Gonzalez et al.l 
(120051) at 2:=0.3,0.5 and 0.9 in the cyan dash-dotted lines, and 
Rus h. Malkan. & Spinogliol (1 19931) at z=0 in the green dash- 
dotted lines. AGNs are excluded from the sample ( il2.2l) . 



iTakeuchi et al.l (l2005h independently estimated the relation 
to be 



log Ltib. = 1.02 + 0.972 log Li2p„ 



(9) 



log Ltir = log(0.89i^;^?) + 1.094 log Li2^^ 



(8) 



which we also use to check our conversion. As both au- 
thors state, these conversions contain an error of factor of 2-3. 
Therefore, we should avoid conclusions that could be affected 
by such eiTors. 

Then the 12/im LF is weighted by the Ltir and integrated 
to obtain TIR density. EiTors are estimated by varying the fit 
within Icr of uncertainty in LFs, and errors in converting from 
Li2^m to Ltir are added. The latter is by far the largest source 
of uncertainty. Best fit parameters are presented in Table |2] In 
FigHni we show total luminosity density based on the 12/im LF 
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Fig. 9. Evolution of 12/j,m IR luminosity density computed by 
integrating the 12/im LPs in Fig|8] 
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Fig. 10. TIR luminosity density computed by integrating the 
12/im LPs in PiglH] 



presented in Pig|8] The resul ts show a rapid increa se of ftjR, 
agreeing with previous work dLe Floc'h et al.L 120051) within the 
errors. 

We also integrate monochromatic Li2fmi over the LPs 
(without converting to Ltir) to derive the evolution of to- 
tal 12/im monochromatic luminosity density, Vli2^m- The re- 
sults are shown in Pig|9] which shows a strong evolution of 



oc (1 + z) 



1.4±1.4 



. It is interesting to compare this to 



oc (1 + 2;)^-^*° '' obtained in il3.2l Although errors are 
significant on both estimates, ^i2^im and fis/im show a possibly 
different evolution, suggesting that the cosmic infrared spectrum 
changes its SED shape. Whether this is due to evolution in dust, 
or dusty AGN contribution is an interesting subject for future 
work. 
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Fig. 11. An example of the SED fit. The red dashed line shows 
the best-fit SED for the UV-optical-NIR SED, mainly to estimate 
photometric redshift. The blue solid line shows the best-fit model 
for the IR SED at A > 6/im, to estimate Ltir- 



3.5. TIRLF 

AKARI's continuous mid-IR coverage is also superior for SED- 
fitting to estimate Ltir, since for star-forming galaxies, the 
mid-IR part of the IR SED is dominated by the PAH emissions 
which reflect the SPR of galaxies, and thus, correlates well with 
Ltir, which is also a good indicator of the galaxy SPR. The 
AKARI's continuous MIR coverage helps us to estimate Ltir- 

After photometric redshifts are estimated using the UV- 
optical-NIR photometry, we fix the redshift at the photo-z, then 
use the same LePhare code to fit the infrared part of the SED 
to estimate TIR luminosity. We used iLagache. Dole. & Pugetl 
(i2003) 's SED templates to fit the photometry using the AKARI 
bands at >6/im (S7, S9W, Sll, L15, L18W and L24). We 
show an example of the SED fit in Pig.IH] where the red dashed 
and blue solid lines show the best-fit SEDs for the UV-optical- 
NIR and IR SED at A > 6/im, respectively. The obtained total 
infrared luminosity (Ltir) is shown as a function of redshift in 
FigEI with spectroscopic galaxies in large triangles. The figure 
shows that the AKARI can detect LIRGs (Ltir > 10"^©) 
up to z=l and ULIRGs (Ltir > lO^^ L©) to z=2. We also 
checked that using di fferent SED models (ICharv & Elbazil200U 
iDale&HelouL 12002 !) does not change our essential results. 

Galaxies in the targeted redshift range are best sampled in 
the 18/im band due to the wide bandpass of the L18W filter 
(iMa tsuhara et al., 2006). In fact, in a single-band detection, the 
18/im image returns the largest number of sources. Therefore, 
we applied the l/T^max method using the detection limit at 
L18W. We also checked that using the LI 5 flux limit does 
not ch ange our main results. The same iLagache. Dole. & Pugei 
(I2003h 's models are also used for fc-corrections necessary to 
compute Knax and Vinin- The redshift bins used are 0.2< 
z <0.5, 0.5< z <0.8, 0.8< z <1.2, and 1.2< z <1.6. A distri- 
bution of Ltir in each redshift bin is shown in PigfT3] 

The obtained Ltir LPs are shown in PiglT4] The uncertain- 
ties are esimated through the Monte Carlo simula t ions ( j l2.4l i. 
Por a local benchmark, we overplot ISanders et alj (l2003h who 
derived LPs from the analytical fit to the IRAS Revised Bright 
Galaxy Sample, i.e., oc L^"-^ for L < L* and oc i^^-^ for 
L > L* with L* = 10^°-^Lq. The TIR LPs show a strong evo- 
lution compared to local LPs. At 0.25 < z < 1.3, L^jj^ evolves 
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Fig. 12. TIR luminosity is shown as a function of photometric 
redshift. The photo-z is estimated using UV-optical-NIR pho- 
tometry. Ltib^ is obtained through SED fit in 7-24/im. 
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Fig. 13. A histogram of TIR luminosity. From low -redshift, 144, 
192, 394, and 222 galaxies are in 0.2< z <0.5, 0.5< z <0.8, 
0.8< z <1.2, and 1.2< z <1.6, respectively. 



Fig. 14. TIR LFs. Vertical lines show the luminosity correspond- 
ing to the flux limit at the central redshift in each redshift bin. 
AGNs are excluded from the sample ( m.2i . 
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as oc (1 + z)"'-^*°-^. We further compare LFs to the previous 
work in Sj4] 

3.6. Bolometric IR luminosity density based on tlie TIR LF 

Using the same methodology as in previous sections, we inte- 
grate Ltir LFs in Fig [14] through a double-power law fit (eq. 
|5] and |6]l. The resulting evolution of the TIR density is shown 
with red diamonds in Fig(T5] which in in good agreement with 
iLe Floc'h et al.l (l2005h within the errors. Errors are estimated by 
varying the fit within Icr of uncertainty in LFs. For uncertainty 
in the SED fit, we added 0. 15 dex of error. Best fit parameters are 
presented in Table |2] In FiglTS] we also show the contributions 
to r^T/i? from LIRGs and ULIRGs with the blue squares and or- 
ange triangles, respectively. We further discuss the evolution of 
nriR in g] 



Fig. 15. TIR luminosity density (red diamonds) computed by 
integrating the total LF in Fig[T4l The blue squares and orange 
triangles are for LIRG and ULIRGs only. 



4. Discussion 

4. 1. Comparison with previous work 

In this section, we compare our results to previous work, espe- 
cially those based on the Spitzer data. Comparisons are best done 
in the same wavelengths, since the conversion from either Lg^m 
or Li2tim to Ltir involves the largest uncertainty. Hubble pa- 
rameters in the previous work are converted to h = 0.7 for com- 
parison. 
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4.1.1. 8/im LFs 

Recently, using the Spitzer space telescope, restfra me E/im LFs 
ofz_~l galaxies have been computed in detail by Can)uti et al.' 
(I2007h in the GOODS fields and by Babbedge et al. (2006) in 
the SWIRE field. In this section, we compare our restframe 8/im 
LFs (FiglDi to these and d iscuss possib l e diffe rences. 

In FigH we overplot ICaputi et al.1 (l2007l) 's LFs at z=l and 
z=2 in the cyan dash-dotted lines. Their z=2 LF is in good agree- 
ment with our LF at 1.8< z <2.2. However, their z=l LF is 
larger than ours by a factor of 3-5 at log£ > 11.2. Note that 
the brightest ends (logL ^ 11.4) are consistent with each other 
to within la. They have excluded AGN using optical-to-X-ray 
flux ratio, and we also have excluded AGN through the optical 
SED fit. Therefore, especially at the faint-end, the contamina- 
tion f rom AGN is not like ly to be the main cause of differences. 
Since ICaputi et al.l (l2007h uses GOODS fields, cosmic variance 
may play a role here. The exact reason for the difference is un- 
known, but we point out that their fljj^ estimate at z=l is also 
higher than other estimates by a factor of a few (see their Fig. 15). 
Once converted int o Ltir, iMagnelli et alj (l2009h also reported 
ICaputi et alJ (l2007h 's z=\ LF is higher than their estimate based 
on 70/zm by a factor of several (see their Fig. 12). They con- 
cluded the difference is from di fferent SED models used, since 
their LF matched with that of iCaputi et alJ (l2007l) 's once the 
same SED models were used. We will compare our total LFs 
to t hose in the litera t ure be low. 

iBabbedge et al.l (l2006h also computed restframe 8/im LFs 
using the Spitzer/SWIRE data. We overplot their results at 
0.25 < z < 0.5 and 0.5 < z < 1 in Fig|4]with the pink dot- 
dashed lines. In both redshift ranges, good agreement is found at 
higher luminosity bins (is^im > 10^° '^L(n). Ho wever, at all red- 
shift r anges including the ones not shown here, iBabbedge et alJ 
(I2006h tends to show a flatter faint-end tail than ours, and a 
smaller by a factor of '^3. Although the exact reason is un- 
known, the deviation starts toward the fainter end, where both 
works approach the flux limits of the surveys. Therefore, possi- 
bly incomplete sampling may be one of the reasons. It is also re- 
ported that the faint-end of IR LFs depends on the environment, 
in the sense that higher-density environment has steeper faint- 
end ta il (Goto et al., 2010). Note that at z=l, Babbedge et al. 
(l2006h 's LF (pink) deviates from that by ICaputi et al.1 (l2007h 
(cyan) by almost a magnitude. Our 8/im LFs are between these 
works. 

These comparisons suggest that even with the current gen- 
eration of satellites and state-of-the-art SED models, factor-of- 
several uncertainties still remain in estimating the 8/zm LFs 
at z^l. More accurate determination has to await a larger 
and deeper survey by the next generation IR satellites such as 
Herschel and WISE. 

To summarise, o ur 8/im LFs are be tween those by 
iBabbedee et"an (l2006l) and lCaputi et al.l (l2007h . and discrepancy 
is by a factor of several at most. We note that both of the previ- 
ous works had to rely on SED models to estimate Lg^m from 
the Spitzer S24fim in the MIR wavelengths where SED model- 
ing is difficult. Here, AKARJ's mid-lR bands are advantageous 
in directly observing redshifted restframe 8/im flux in one of the 
AKARl's filters, leading to more reliable measurement of 8/im 
LFs without uncertainty from the SED modeling. 

4.1.2. 12/im LFs 

iPerez-Gonzalez et alj 42005") investigated the evolution of rest- 
frame 12/im LFs using the Spitzer CDF-S and HDF-N data. We each galaxy. 



overplot their results in similar redshift ranges as the cyan dot- 
dashed lines in FiglH] Considering both LFs have significant er- 
ror bars, these LFs are in good agreement with our LFs, and 
show significant evolution in the 12/zm LFs compa red with the 
z=0 12^mLFbv lRush. Malkan. & Spinogliol(ll993h . The agree- 
ment is in a stark contrast to the comparison in 8/im LFs in 
fjJTTT] where we suffered from differnces of a factor of several. 
A possible reason for this is that 12/im is sufficiently redder than 
8/im, that it is easier to be extrapolated from S24^m in case of 
the Spitzer work. In fact, at z=\, both the Spitzer 24/im band 
and AKARI L24 observe the restframe 12/im directly. In addi- 
ton, mid-IR SEDs around 12/im are flatter than at 8/im, where 
Y'AH emissions are prominent. Therefore, SED models can pre- 
dict the flux more accurately . In fa ct, this is part of the rea- 
son why'Perez-Gonzalez et al.'('2005') c hose to in vestigate 12/im 
LFs. Perez-Gonzalez et al. (2005) used lCharv & Elbaz (2001)'s 
SED to extrapolate 524pm, and yet, they agree well with AKARI 
results, which are derived from filters that cover the restframe 
12/im. However, in other words, the discrepancy in 8/im LFs 
highlights the fact that the SED models are perhaps still imper- 
fect in the 8/im wavelength range, and thus, MIR-spectroscopic 
data that covers wider luminosity and redshift ranges will be 
necessary to refine SED models in the mid-I R. AKARl's mid- 
IR slitless spectroscopy survey dWadal l2008h may help in this 
regard. 



4.1.3. TIR LFs 

Lastly, we compare our TIR LFs (FigfT4l) with those in the liter- 
ature. Although the TIR LFs can also be obtained by converting 
8/im LFs or 12/i,m LFs, we already compared our results in these 
wavel engths in the last subse ctio ns. Here, we compare our TIR 
LFs to lLe Floc'h et alj (l2005h and lMagneUi et al.1 (l2009l) . 

i Le Floc'h et~an (l2005h obtained TIR LFs using the Spitzer 
CDF-S data. They have used the best-fit SED among various 
templates to estimate Ltir.- We overplot their total LFs in FigfT4l 
with the cyan dash-dotted lines. Only LFs that overlap with our 
redshit ranges are shown. The agreement at 0.3 < z < 0.45 
and 0.6 < z < 0.8 is reasonable, considering the error bars on 
both s ides. However, in all three redshift ranges, Le Floc'h et al] 
( 120051) 's LFs are higher than ours, especially for 1.0 < z < 1.2. 

We also overplot TIR LFs by 'Magnelli et al.' ('2009'), who 
used Spitzer 70/im flux and Charv & Elbaz (2001)'s model to 
estimate Ltir- In the two bins (centered on 2=0.55 and z=0.85; 
pink dash-dotted lines) which closely overlap wit h our redshift 
bins, e xcellent agreement is found. We also plot iHuvnh et al.l 
('200V)'s LF at 0.6 < z < 0.9 in the navy dash-dotted fines, 
which is computed from Spitzer 70/im imaging in the GOODS- 
N, and this also shows very good agreement with ours. These 
LFs are on top of each other within the eiTor bars, despite the 
fact that these measurements are from different data sets using 
different analyses. 

This m eans that Le Floc'h et alj (l2005h 's LFs is also higher 
than that o f lMagnelli et alJ(l2009l) in add ition to ours. A possible 
reason is that both lMagnelli et al.' ("2009*) and we removed AGN 
(at least bright ones), whereas Le Floc'h et al. (2005) included 
them. This also is consistent with the fact that the difference 
is larger at 1.0 < z < 1.2 where both surveys are only sen- 
sitive to luminous IR galaxies, which are domi nated by AGN 



Anoth er possible source of uncertainty is that Magnelli et al 



(2009) and we used a single SED library, while Le Floc'h et al. 



20051) picked the best SED template among several libraries for 
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Fig. 16. Evolution of TIR luminosity density based on TIR LPs (red circles), 8/im LPs (stars), and 12/im LPs (filled triangles). The 
blue open squares and orange filled sq uares are for LIRG and ULIRGs only, also bas ed on our Lttr LPs. Overplott e d dot-dashed 
lines are estimates from the literature: iLe Floc'h et alj (l2005h . iMagnelh etal] (|2009|) , iPerez-Gonzdlez et all (l2005l) . ICaputi et all 
(l2007h . andlBabbedae et al. (200j ) are in cyan, yellow, green, navy, and pink, respectively. The purple dash-d otted line shows UV 
estimate bv lSchiminovich et al.l (l2005h . The pink dashed line shows the total estimate of IR (TIR LP) and UV (ISchiminovich et all 
I2OOI . 



4.2. Evolution ofiliR 

In this section, we compare the evolution of fi/ ^ as a function 
of redshift. In Pig[T6] we plot ^Im estimated from TIR LPs (red 
circles), 8/im LPs (brown stars), and 12/im LPs (pink filled tri- 
angles), as a function of redshift. Estimates based on 12/im LPs 
and TIR LPs agree each other very well, while those from 8//m 
LPs show a slightly higher value by a factor of a few than oth- 
ers. This perhaps reflects the fact that S/zm is a more difficult 
part of the SED to be modeled, as we had a poorer agreement 
among papers in the literature in 8/im LPs . The bright-end slope 
of the double-power law was 3.5+°- ^ in Tabled This is flat- 



ter than a Schecht er fit by[B abbedge et alj (120061) and a double- 
exponential fit bv ICaputi et alJ (l2007ij) . This is perhaps why we 
obtained higher ilm in 8/xm. 

We overplot estimates from various papers in the litera- 
ture jLe P loc'h et al.l '2005'; 'Babbedge et al.', '2006'; 'Capu ti et all 
l2007t iPerez-Gonzalez et al., 2005; Magnelh et al., 2009) in the 
dash-dotted lines. Our iljj^ has very good agreement with these 
atO < z < 1.2, with almost all the dash-dotted lines lying within 
our error bars of from Ltib. and 12/im LPs. This is perhaps 
because an estimate of an integrated value such as fi/ ^ is more 
reliable than that of LPs. 

A t z > 1.2, our UiR sho ws a hint of con t inuou s increase, 
while ICaputi et all (l2007[) and lBabbedge etalJ (120061) observed 



a slight decline at z > 1. However, as both authors also pointed 
out, at this high-redshift range, both the AKARI and Spitzer 
satellites are sensitive to only LIRGs and ULIRGs, and thus the 
extrapolation to fainter luminosities assumes the faint-end slope 
of the LPs, which could be uncertain. In addition, this work has 
a poorer photo-z estimate at z > 0.8 (^^^=0.10) due to the rel- 
atively shallow near-IR data. Several authors tried to overcome 
this problem by stacking undetected sources. However, if an un- 
detected source is also not detected at shorter wavelengths where 
positions for stacking are obtained, it would not be included in 
the stacking either. Next generation s atellite such as Herschel, 
WISE, and SPICA (Nakagawal l2008l) will determine the faint- 
end slope at z > 1 more precisely. 

We parameterize the evolution of ftjij using a following 
function. 



niR{z) cx (i + z)^ 



(10) 



By fitting this to the Hir from TIR LPs, we obtained 7 = 
4.4 ± 1.0. T his is consiste n t with most previous works. 
Por example. iLe Floc'h et all (|2005|) obtained 7 = 3.9 ± 
0.4, Perez-Gonzalez et all (l2005l) obtained 7 = 4.0 ± 0.2 
Babbedge et al. (2006) obtained 7 = 4.5l° g, iMagnelh et aP 
(2009) obtained 7 = 3.6 ± 0.4. The agreement was expected 
from PiglT6| but confirms a strong evolution of fi/fl. 
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Fig. 17. Contribution of rtxi r to fttotai 
shown as a function of redshift. 
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4.3. Differential evolution among ULIRG, LIRG, normal 
galaxies 

In Fig. [15] we also plot the contributions to rtjR from LIRGs 
and ULIRGs (measured from TIR LFs) with the blue open 
squares and orange filled squares, respectively. Both LIRGs 
and ULIRGs show strong evolution, as has been seen for to- 
tal iljR in the red filled circles. Normal galaxies {Ltir < 
IO^^Lq) are still dominant, but decrease their contribution to- 
ward higher redshifts. In contrast, ULIRGs continue to increase 
their contribution. From z=0.35 to z=1.4, fljR by LIRGs in- 
creases by a factor of ^1.6, and ftiR by ULIRGs increases by 
a factor of ~10. The physical origi n of ULIRGs in the local 
Universe is often merger/interaction ( Sanders & Mirabeli Il996t 
iTaniguchi & Shioval [l998'; 'Goto. '2005*). It would be interesting 
to investigate whether the merger rate also increases in propor- 
tion to the ULIRG fraction, or if different mechanisms can also 
produce ULIRGs at higher redshift. 



4.4. Comparison to the UV estimate 

We have been emphasizing the importance of IR probes of the 
total SFRD of the Universe. However, the IR estimates do not 
take into account the contribution of the unabsorbed UV light 
produced by the young stars. Therefore, it is important to esti- 
mate how significant this UV c ontribution is. 

[Schiminovich et al. (l2005h found that the energy density 
measured at 1500A evolves as oc (1 + z)'^-5±f^-^ at < z < 1 
and oc (1 + 2:)°-5='=° '' at z > I. using the GALEX data sup- 
plemented by the VVDS spectroscopic redshifts. We overplot 
their UV estimate of psFR. with the purple dot-dashed line in 
Fig(T6] The UV estimate is almost a factor of 10 smaller than 
the IR estimate at most of the redshifts, confirming the impor- 
tance of IR probes when investing the evolution of the total cos- 
mic star formation density. In Fig(T6]we also plot total SFD (or 
^totai) by adding fluv and Uti r, with the magenta dashed line. 
In FigllT] we show the ratio of the IR contribution to the to- 
tal SFRD of the Universe {UtirI ^tir + ^uv) as a function 
of redshift. Although the errors are large, FigllT] agrees with 



iTakeuchL Buat. & Burgarellal (l2005h . and suggests that flriR 
explains 70% of fltotai at z=0.25, and that by z=1.3, 90% of 
the cosmic SFD is explained by the infrared. This implies that 
Utir provides good approximation of the fttotai at z > 1. 

5. Summary 

We have estimated restframe 8/im, llfim, and total infrared lu- 
minosity functions using the AKARI NEP-Deep data. Our ad- 
vantage over previous work is AKARI's continuous filter cov- 
erage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, 
and 24/im), which allow us to estimate mid-IR luminosity with- 
out a large extrapolation based on SED models, which were the 
largest uncertainty in previous work. Even for Ltir, the SED 
fitting is much more reliable due to this continuous coverage of 
mid-IR filters. 

Our findings are as follows: 

- 8/im LFs show a strong and continuous evolution from 
2=0.35 to z=2.2. Our LFs are la rger than iBabbedge et al.l 
( |2006l ), but smaller thanl Caputi et al. (2007). The difference 
perhaps stems from the different SED models, highlighting 
a difficulty in SED modeling at wavelengths crowded by 
strong R\H emissions. L^^^ shows a continuous evolution 



as L 



OC (1 +z) 



1.6±0.2 ; 



in the range of 0.48 < z < 2. 



- lliim LFs show a strong and continuous evolution from 

z=0.15 to z=1.16with L*^„^ oc (l + z )i-5±o.4^ xhis agrees 

' — I J — — I 

Perez-Gonzalez et al.l (120051) . including a flatter 

faint-end slope. A better agreement than with 8/im LFs was 
obtained, perhaps because of smaller uncertainty in model- 
ing the 12/im SED, and less extrapolation needed in Spitzer 
24/im observations. 

- The T IR LFs show good agr eement with Magnelli et al.l 



(I2009I) . but are smaller than iLe Floc'h et alJ (l2005h . M 



0.25 < z < 1.3, i^^fl evolves as oc (1 + z)'^-'^'^°-^. Possible 
causes of the disagreement include different treatment of 
SED models in estimating Ltir, and AGN contamination. 
- TIR densities estimated from 12/im and TIR LFs show a 
strong evolution as a function of redshift, with ftiR oc 



(1 



+ z 



,4.4±1.0 



^Iir{z) also show an excellent agreement 



with previous work at z < 1.2. 

- We investigated the differential contribution to ftiR by 
ULIRGs and LIRGs. We found that the ULIRG (LIRG) con- 
tribution increases by a factor of 10 (1.8) from z=0.35 to 
z=1.4, suggesting IR galaxies are more dominant source of 
rtiR at higher redshift. 

- We estimated that fli r captures 80% of the cosmic star for- 
mation at redshifts less than 1, and virtually all of it at higher 
redshift. Thus adding the unobscured star formation detected 
at UV wavelengths would not change SFRD estimates sig- 
nificantly. 
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